Perturbative — nonperturbative connection in quantum mechanics and field theory 
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On the occasion of this ArkadyFest, celebrating Arkady Vainshtein's 60 th birthday, I review some 
selected aspects of the connection between perturbative and nonperturbative physics, a subject to 
which Arkady has made many important contributions. I first review this connection in quantum 
mechanics, which was the subject of Arkady's very first paper. Then I discuss this issue in relation to 
effective actions in field theory, which also touches on Arkady's work on operator product expansions. 
Finally, I conclude with a discussion of a special quantum mechanical system, a quasi-exactly solvable 
model with energy-reflection duality, which exhibits an explicit duality between the perturbative and 
nonperturbative sectors, without invoking supersymmetry. 
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DIVERGENCE OF PERTURBATION THEORY 



"The majority of nontrivial theories are seemingly unstable at some 
phase of the coupling constant, which leads to the asymptotic 
nature of the perturbative series." 
A. Vainshtein, 1964 [1] 



In this talk I review some aspects of the historical development of the connection between perturbative and nonper- 
turbative physics. It is particularly appropriate to look back on this subject on the occasion of Arkady Vainshtein's 
60 th birthday, because this has been a central theme of many of Arkady's great contributions to theoretical physics. 
In fact, in his very first physics paper [1], now almost 40 years ago, Arkady made a fundamental contribution to 
this subject. This paper was published as a Novosibirsk report and so has not been widely circulated, especially in 
the West. For this ArkadyFest, Misha Shifman has made an English translation of this paper, and both the original 
Russian and the translation are reprinted in these Proceedings. 

The physical realization of the possibility of the divergence of perturbation theory is usually traced back to a 
profound and influential paper by Dyson [2], in which he argued that QED perturbation theory should be divergent. 
Dyson's argument goes like this: a physical quantity in QED, computed using the standard rules of renormalized 
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QED perturbation theory, is expressed as a perturbative series in powers of the fine structure constant, a = f-: 



F{e 2 ) = c + c 2 e 2 + c 4 e 4 + 



(1) 



Now, suppose that this perturbative expression is convergent. This means that in some small disc-like neighborhood 
of the origin, F(e 2 ) has a well-defined convergent approximation. In particular, this means that within this region, 
F(—e 2 ) also has a well-defined convergent expansion. Dyson then argued on physical grounds that this cannot be 
the case, because if e 2 < the vacuum will be unstable. This, he argued, is because with e 2 < like charges attract 
and it will be energetically favorable for the vacuum to produce e + e~ pairs which coalesce into like-charge blobs, a 
runaway process that leads to an unstable state: 



"Thus every physical state is unstable against the spontaneous creation 
of large numbers of particles. Further, a system once in a pathological 
state will not remain steady; there will be a rapid creation of more 
and more particles, an explosive disintegration of the vacuum 
by spontaneous polarization." 
F. J. Dyson, 1952 [2] 



The standard QED perturbation theory formalism breaks down in such an unstable vacuum, which Dyson argued 
means that F(—e 2 ) cannot be well-defined, and so the original perturbative expansion (1) cannot have been convergent. 

Dyson's argument captures beautifully an essential piece of physics, namely the deep connection between instability 
and the divergence of perturbation theory. The argument is not mathematically rigorous, and does not prove one way 
or another whether QED perturbation theory is convergent or divergent, or analytic or nonanalytic. However, it is 
nevertheless very suggestive, and has motivated many subsequent studies in both quantum mechanics and quantum 
field theory. 

At roughly the same time, C. A. Hurst [3] and W. Thirring [4] (see also A. Petermann [5]) showed by explicit 
computation that perturbation theory diverges in scalar <j) 3 theory. Both Hurst and Thirring found lower bounds on 
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the contribution of Feynman graphs at a given order of perturbation theory, and showed that these lower bounds were 
themselves factorially divergent. Hurst used the parametric representation of an irreducible, renormalizcd and finite 
<f> 3 Feynman graph, to show that the magnitude of this graph was bounded below: 



l/l > 



n -n+3/2-E/2 £ 2n-2 n -(n-E+5)/2 
{2ir) F 2™+!/2-3iJ/2 3 4n-7/2 7 (n+B-4)/2 



A™ 



(2) 



Here n is the loop order, A is the cubic coupling constant, E is the number of external lines, F = \{Sn — E) is the 
number of internal lines, and 7 is a constant depending on the external momenta. This lower bound is found by clever 
rearrangements of the parametric representation, together with the identity 
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The second important piece of the argument is to show that there are no sign cancellations which would prevent this 
lower bound from a typical graph from being used to obtain a lower bound on the total contribution at a given order. 
This requires some technical caveats - for example, for a two-point function one requires p 2 < to 2 . The final piece of 
Hurst's argument is the fact that the number of distinct Feynman diagrams at n th loop order grows like (^)\n\. 

Together, the lower bound (2), the nonalternation of the sign, and the rapid growth of the number of graphs, lead 
to a lower bound for the total contribution at n-loop order (with E external lines): 



^2 1 > C" n n/2+5/2 ~ E/2 A" 



(4) 



Here C is a finite constant, independent of n. Therefore, Hurst concluded that in A<p 3 theory, perturbation theory 
diverges for any coupling A. He also suggests that a similar argument should hold for A0 4 theory, and comments: 



"If it be granted that the perturbation expansion does not lead to a 
convergent series in the coupling constant for all theories which can 
be renormalized, at least, then a reconciliation is needed between this 
and the excellent agreement found in electrodynamics between 
experimental results and low-order calculations. It is suggested that this 

agreement is due to the fact that the S-matrix expansion is to be 
interpreted as an asymptotic expansion in the fine-structure constant ..." 

C. A. Hurst, 1952 [3] 



Thirring's argument [4] was similar in spirit, although he concentrated on the <fi 3 self-energy diagram. Thirring 
found a set of graphs that were simple enough that their contribution could be estimated and bounded below, while 
plentiful enough that they made a divergent contribution to the perturbative series. He noted that the proof relied 
essentially on the fact that certain terms always had the same sign, and traced this fact to the hermiticity of the 
interaction. He found the following (weaker) lower bound, valid for p 2 
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Thirring concluded that there was no convergence for any A. His final conclusion was rather pessimistic: 



"To sum up, one can say that the chances for quantized fields to 
become a mathematically consistent theory are rather slender." 
W. Thirring, 1953 [41 



These results of Dyson, Hurst and Thirring, provide the backdrop for Arkady's first paper [1], "Decaying systems and 
divergence of perturbation theory" , written as a young student beginning his PhD at Novosibirsk. I encourage the 
reader to read Arkady's paper - it is simple but deep. I paraphrase the argument here. The main contribution of his 
paper was to provide a quantitative statement of the relation between the divergence of perturbation theory and the 
unstable nature of the ground state in (j) 3 theory. 

Motivated by the earlier results for (j) 3 theory (in 4 dimensions), Arkady had the clever idea to consider </> 3 theory 
in + 1 dimensions, which is just quantum mechanics. Here it is natural to consider the hamiltonian 



H=\4> 2 + \m 2 cf-\4? 



(6) 



FIG. 1: Unstable ground state for V{(j>) = \<t> 2 - \<f A 



and the ground state \ip) such that H\ip) = E\ip). To make connection with the field theory results, note that the 
two-point function 



iG(t) 



is related to the energy E as 
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where in the last step we have used the fact that, by dimensional reasoning, the energy E can be expressed as 
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E = mf( — 5). Thus, if the pcrturbativc expression for the two-point function diverges, the expression for the ground 
state energy, E = E(X 2 ), should also diverge. One subtlety here is that the state \ip) is clearly unstable. Arkady 
showed in an appendix [1] how to deal with this, by considering the adiabatic evolution of a stable state into an 
unstable state. In particular this suggests that the expression for E = E(\ 2 ) must have a cut along the positive A 2 
axis, as shown in Fig. 2, with an associated jump in the imaginary part of E across this cut. 




FIG. 2: The complex z = A 2 plane, showing the cut along the positive 2-axis. 

Under the (important) assumption that there are no other cuts or poles in the complex z = A 2 plane, Cauchy's 
theorem implies that : 
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Thus, the perturbative expansion coefficients are explicitly related to the moments of the imaginary part of the 
energy along the cut. Furthermore, it is clear from (9) that at large n (i.e., at large order in perturbation theory), the 
dominant contribution comes from the behavior of ImE(z) as z — > 0. This observation is very important, because 
the z = A 2 — > limit is a semiclassical limit (note that the barrier height goes like j?, and the barrier width like j). 
Hence, in this limit the imaginary part of the energy may be estimated using semiclassical techniques, such as WKB. 
Simple scaling shows that the leading WKB approximation for the imaginary part of the energy has the form 

ImE(z) ~ -^=e- h / z , z^O (10) 

where a and b > are (calculable) constants. Note, of course, that this expression is nonperturbative in z = A 2 . 
Now consider the perturbative expansion for the lowest energy 

E(X 2 )= 1 - + J2 c ^ 2n (11) 

Z n=l 

Inserting the WKB estimate (10) into the dispersion relation (9), we see that at large order the perturbation theory 
coefficients should behave as 
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So, this argument suggests that the perturbative expansion (11) for E(\ 2 ) should be a divergent nonalternating series. 
Indeed, it is straightforward to do this perturbative calculation to very high orders, and to do the WKB calculation 
precisely, and one finds excellent agreement [6]. If the hamiltonian is rescaled as H = p 2 + jx 2 — Xx 3 (this scaling 
makes the expansion coefficients integers), then the leading growth rate for large n is 
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(13) 



which agrees with Arkady's form (12), and fits beautifully the growth rate of the actual expansion coefficients [6]. 
Indeed, the factorial growth of the perturbative coefficients kicks in rather early, as illustrated in Fig. 3. 
The most important physics lessons from Arkady's paper [1] are : 

(i) the divergence of perturbation theory is related to the possible instability of the theory, at some phase of the 
coupling. 

(ii) there is a precise quantitative relation (9) between the large-order divergence of the perturbative coefficients 
and nonperturbative physics. 
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FIG. 3: The ratio of the exact perturbative coefficients c n in (11) to the leading WKB expression from (13), as a function of 
the order n. 

In more modern language, this divergence associated with instability and tunneling is a divergence due to instantons. 
This idea has become a cornerstone of quantum field theory [7]. However, since this time, it has been found that in 
quantum field theory (as distinct from quantum mechanics) there is yet another source of divergence in perturbation 
theory - this divergence is due to " renormalons" , which arise essentially because of the running of the coupling 
constant, and can be related to special classes of diagrams [8]. For recent developments in this important subject, see 
the talks by M. Benekc, I. Balitsky and E. Gardi in these Proceedings. 

The connection between the large-order behavior of quantum mechanical perturbation theory and WKB methods 
was developed independently, and was probed in great depth, by Bender and Wu, who studied the quartic anharmonic 
oscillator [9] . Bender and Wu developed recursion techniques for efficiently generating very high orders of perturbation 
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FIG. 4: Stable ground state for the potential V = x 2 + \x 4 , on the left, and unstable ground state for the potential V = x 2 — Xx 4 
on the right. 



theory, and compared these results with higher orders of the WKB approximation. Their WKB analysis is a tour de 
force, and the agreement with the large order perturbative coefficients is spectacular. 

The dependence of the ground state energy on the coupling A in the X(j) 4 case is different from the A</> 3 case. As is 
clear from Fig. 4, in the \<fi case the instability arises when the coupling A changes sign from positive to negative. 
Thus, the cut is expected along the negative real axis in the complex z = A plane, which led Bender and Wu to the 
following dispersion relation: 
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(14) 



where F(X) = j(E(X) — i) has one subtraction. Here, Bender and Wu made use of some rigorous results [10] 
concerning the analyticity behavior of E(X) in the complex A plane : (i) |£7(A)| <~ |A| 1//3 for large |A|, (ii) E(X) is 
analytic in the cut A-plane, with the cut along the negative real axis, and (hi) the expansion E(X) <~ \ + J2^=i c nX n 
is asymptotic. 

The dispersion relation (14) relates the perturbative expansion coefficients for F(X) to the discontinuity of F(X) 
across the cut. Bender and Wu used high orders of WKB to compute this nonperturbative imaginary part, thereby 
providing an extremely precise connection between the large orders of perturbation theory and semiclassical tunneling 
processes. They found that the expansion coefficients are alternating and factorially divergent 
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Bender and Wu extended these results to the general anharmonic oscillator: 

-^ + ^ + A^_ £(A) ), = o 

Then the k th energy level has an asymptotic series expansion 

1 n=l 



(15) 
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where the perturbative expansion coefficients c^ N are related to the lifetime of the k th unstable level when the coupling 
A is negative: 
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(18) 



Once again, these growth estimates are derived from WKB and fit the actual perturbative coefficients with great 
precision [9]. 

An important distinction between the A^> 3 and A</> 4 quantum mechanical oscillators is that the perturbative series 
for the energy eigenvalue is nonalternating in the A</> 3 case, and alternating in the X<p 4 case. This is directly related 
to the fact that the A</> 3 case is inherently unstable (for any real A), while the A0 4 case is stable if A > 0, but unstable 
if A < 0. Even though the quartic anharmonic oscillator with A > is completely stable, the perturbative expression 
for the ground state energy is divergent, because the theory with A < is unstable. 
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It is interesting to note that there has been much recent interest in certain nonhermitean hamiltonians, whose 
spectra appear to be completely real, despite the nonhcrmicity [11]. For example, there is very strong numerical 
evidence [12, 13] that the hamiltonian 



H =p 2 



1 



(19) 



has a completely real spectrum. For the massless case, H = p 2 



the reality of the spectrum has in fact been 



proved rigorously [14]. These results are nicely consistent with Arkady's analysis of the \<fr 3 quantum mechanical 
model, since an imaginary coupling corresponds to z = A 2 being real and negative. In this region it was assumed that 
E{z) is analytic (recall from Fig. 2 that the cut is along the positive z axis), and the perturbative series is divergent 
but alternating, and can be analyzed using various standard (Pade and Borcl) techniques [12, 13], yielding excellent 
agreement with numerical integration results. 

A useful mathematical technique for dealing with divergent series is Borel summation [10, 15-17]. This method is 
best illustrated by the paradigmatic case of the factorially divergent series. The series 



/( 5 ) ~ ]T(-i)W 



(20) 



is clearly divergent, and is alternating if g > 0. Using the standard integral representation, n! = dse s s n , and 
formally interchanging the summation and integration, we can write 
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ds ■ 



-s/g 



1 + 8 



(21) 



This integral representation is defined to be the Borel sum of the divergent series in (20). The advantage of the 
integral is that it is convergent for all g > 0. To be more precise, all this actually shows is that the integral in (21) 
has the same asymptotic series expansion as the divergent series in (20). In order for this identification between the 
series and the Borel integral to be unique, various further conditions must be satisfied [15-17]. In some quantum 
mechanical examples it is possible to study these conditions rigorously [10], but unfortunately this is usually impossible 
in realistic quantum field theories. This means we arc often confined to "experimental mathematics" when applying 
Borel techniques to perturbation theory in QFT. Nevertheless, I prefer the attitude of Heaviside: 



"The series is divergent; therefore we may be 
able to do something with it" 
O. Heaviside, 1850 - 1925 



to the (older) attitude of Abel: 



"Divergent series are the invention of the devil, and it is shameful to 
base on them any demonstration whatsoever" 
N. H. Abel, 1828 



Continuing the paradigm in (20), when g < the series (20) becomes nonalternating. Then the same formal steps 
lead to the following representation: 
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ds ■ 



-s/g 



(22) 



Clearly there is a problem here, as there is a pole on the contour of integration, and so an ambiguity enters in the way 
one treats this pole. This means that the nonalternating factorially divergent series in (20), with g < 0, is not Borel 
summable. For example, the principal parts prescription leads to the following imaginary part for the nonalternating 
series : 



Imf(-g) = - exp 



(23) 



This imaginary contribution is nonperturbative in the expansion parameter g. It is not seen at any finite order 
in perturbation theory. However, the imaginary part is inherently ambiguous in the absence of further physical 
information beyond the series expansion (20) itself. 
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Despite this ambiguity, it should be clear that the Borel approach provides a natural formalism with which to analyze 
the problem of the divergence of perturbation theory It captures the essence of the connection with nonpcrturbative 
tunneling, and associates such nonperturbative effects with the unstable cases of the A</> 3 oscillator and the \<fi 4 
oscillator with A < 0, for which the perturbative series is indeed nonalternating and factorially divergent. 

Similar formal expressions exist for the Borel sum of f(g) ~ c n9 n , if the c n are not simply factorial as in (20), 
but have the general form: 



c„~ (3 n Y( in + 5) 

where (3, 7 > and 8 are constants. Then the Borel sum approximation is 
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For g < the nonalternating series has an imaginary part: 
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In the next section we will use these relations in an explicit example. 

Much more could be said about the divergence of perturbation theory, both in quantum mechanics and field theory. 
Lipatov [18] generalized the instanton technique to scalar field theory, showing that large orders of perturbation theory 
may be described by pseudoclassical solutions of the classical field equations, together with quantum fluctuations. 
This approach built on instanton results of Langer [19] in his classic study of metastability. Perturbation theory 
for systems with degenerate minima was investigated in [20], and a new twist on the perturbative-nonperturbative 
connection for this degenerate case is discussed in section 3 of this talk. For further references, I refer the interested 
reader to the review [21] of Le Guillou and Zinn-Justin as an excellent source. 



II. EFFECTIVE ACTIONS, OPES, AND DIVERGENT SERIES 

As mentioned previously it is extremely difficult, even in quantum mechanics, to prove truly rigorous results 
concerning the divergence of perturbation theory [10]; in quantum field theory we are even more restricted when 
it comes to rigor. However, the study of effective actions is an example in QFT where some rigorous results are 
possible. This also makes connection with the subject of operator product expansions, which is another subject to 
which Arkady has made seminal contributions. 

For this talk, I consider the QED effective action, which encodes nonlinear interactions due to quantum vacuum 
polarization effects, such as light-by-light scattering. The effective action is defined via the determinant that is 
obtained when the electron fields are integrated out of the QED functional integral: 

S[A] =-^logdet(^ 2 + m 2 ) (27) 

where Ip = "f^D^ = 7 M (9 M — ieA^) is the Dirac operator in the classical gauge field background A^. The effective 
action has a natural perturbative expansion in terms of the electromagnetic coupling e, as represented in Fig. 5. 
Indeed, by charge conjugation invariance (Furry's theorem), the expansion involves only even numbers of external 




FIG. 5: Perturbative expansion of the one-loop effective action 

2 

photon lines, which means that the perturbative series is actually a series in the fine structure constant a = 
Another natural expansion is the "effective field theory" expansion (or "large mass" expansion): 

S[A]=m*'£c n ¥— (28) 

n 
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Here represents gauge invariant and Lorentz invariant terms constructed from the field strength F^ v , and having 
mass dimension n. For example, at mass dimension 8, we can have (F^F^) 2 or (F^F^) 2 , while at mass dimension 
10 we could have a term {d^F vp d ti F vp ){F a pF OL ^). As shown by Arkady and his collaborators [22, 23], such an 
expansion is related to the operator product expansion (OPE), such as that for 



(29) 



In the special case where the classical background has constant field strength F^, the perturbative and large mass 
expansions coincide. This case of a constant background field was solved by Eulcr and Heisenberg [24] (see also 
Weisskopf [25] and Schwinger [26]), who obtained an exact nonperturbative expression for the effective action: 



S[A] = ii"$ e ~"" 2s { (es)2|g|cot es (V^ 2 +£ 2 +^) 



x coth 



(V^ 2 + Q 2 - t) 



(30) 



Here T = \F liV F lLV , and Q — jF^F^ u , are the two Lorentz invariant combinations. The —1 term in the integrand 
corresponds to the zero- field subtraction of S[0], while the last term, |(es) 2 .F, corresponds to charge renormalization 
[26]. There are several important physical consequences of this result [24-26]. First, expanding to leading order in 
the fields, we find the famous light-by-light term: 
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(31) 



Second, for a constant electric field background, there is an imaginary part 
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(32) 



which gives the pair production rate due to vacuum polarization. 

The Euler-Heisenberg result (30) provides an excellent example of the application of Borel summation techniques, 
as we now review. Consider first of all the case of a uniform magnetic field background of strength B. Then the full 
perturbative expansion of the Euler-Heisenberg result is 
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(33) 



Viewed as a low energy effective action, the "low energy" condition here is simply that the characteristic energy 
scale for electrons in the magnetic background, K—, is much smaller than the electron rest mass scale mc 2 . The 
expansion coefficients in the series (33) involve the Bernoulli numbers B2n, which alternate in sign and grow factorially 
in magnitude [27]. Thus, the series (33) is an alternating divergent series. In fact, the expansion coefficients are: 
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(34) 



If we keep just the leading term in (34), then the expansion coefficients are of the form in (24), so that the Borel 
prescription (25) yields the leading Borel approximation for S as 
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(35) 



In Fig. 6 this leading Borel approximation is compared with successive terms from the perturbative expansion (33). 
Clearly the Borel representation is far superior for larger values of the expansion parameter ^4. 

Actually, in this Euler-Heisenberg case, we can do even better since we are in the unusual situation of having the 
exact expression (34) for the perturbative coefficients, to all orders. Furthermore, the subleading terms in (34) are 
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FIG. 6: Comparison of the leading Borel expression (35) [short-long-dash curve], as a function of with the exact expression 
(36) [solid curve], and successive partial sums from the series (33) [short-dash curves]. The leading Borel expression is much 
better than the series expressions for ^4- > 1. 
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also of the form (24), for which we can once again use the Borel prescription (25). Including these subleading terms, 
we find 



It is interesting to note that the correct renormalization subtractions appear here. The integral representation (36) is 
precisely the Euler-Heisenberg expression (30), specialized to the case of a purely magnetic field. Thus, the "proper- 
time" integral representation (36) is the Borel sum of the divergent perturbation series (33). Conversely, the divergent 
series (33) is the asymptotic expansion of the nonperturbative Euler-Heisenberg result (30). 

Now consider the case of a purely electric constant background, of strength E. Pcrturbatively, the only difference 
from the constant magnetic case is that we replace B 2 — ► — E 2 . This is because the only Lorentz invariant combination 
is (E 2 — B 2 ) (clearly, E ■ B = if one or other of E or B is zero). Thus, the alternating series in (33) becomes 
nonalternating, without changing the magnitude of the expansion coefficients. Applying the Borel dispersion relation 
(26) leads to an imaginary part in exact agreement with the nonperturbative pair-production result (32). 

The divergence of the Euler-Heisenberg effective action is very well known [28-32]. This example is analogous to 
Dyson's argument. The Euler-Heisenberg perturbative series cannot be convergent, because if it were convergent, then 
for very weak fields there would be no essential difference between the magnetic or electric background. However, we 
know from nonperturbative physics that there is an inherent vacuum instability in the case of an electric background, 
and that this leads to an exponentially small imaginary part to the effective action, which corresponds to pair 
production due to vacuum polarization. 

The phenomenon of pair production from vacuum is fascinating, but is very difficult to observe because the expo- 
nential suppression factor is exceedingly small for realistic electric field strengths. The critical electric field at which 
the exponent becomes of order 1 is E c = ~ 10 16 V/cm, which is still several orders of magnitude beyond the 

peak fields obtained in the most intense lasers. See [33] for a recent discussion of the prospects for observing vacuum 
pair production using an X-ray free-electron laser. 

Even if one can produce such an intense background electric field, it is clear that it will not be a uniform field. Thus, 
it is important to ask how the Eulcr-Hciscnberg analysis is modified when the background field is inhomogeneous. 
This is a very difficult problem in general. The most powerful approach is through semiclassical WKB techniques 
[34, 35]. However, here I turn this question around and ask what this issue can tell us about the series expansion 
of the effective action when the background is inhomogeneous [36]. If the field strength is inhomogeneous, then the 
large mass expansion (28) of the effective action involves many more terms at a given order, since we can now include 
terms involving derivatives of F^. In fact [37], the number of terms appears to grow factorially fast: 1, 2, 7, 36, ... 
. But there are two obvious problems with quantifying this divergence. First, the expansion is not a true series, for 
the simple reason that at successive orders many completely new types of terms appear. Second, to learn anything 
nonperturbative one needs to go to very high orders in the derivative expansion, which is extremely difficult. 




(36) 



where we have used the trigonometric expansion [27] : 




(37) 
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Fortunately, there is a special case in which both these problems are circumvented at once [36]. Consider the 
particular inhomogcncous magnetic field pointing in the z-dircction, but varying in the rr-dircction as: 



B(x) = Ssech 2 (x/A) 



(38) 



where A is an arbitrary inhomogeneity scale. In such a background the effective action can be computed in closed- 
form using nonperturbative techniques [38]. Also, since the inhomogeneity of the background is encoded solely in the 
dependence on the scale A, the large mass expansion (28) can be expressed as a true series. That is, at a given order 
in the derivative expansion, all terms with the same number d of derivatives combine to give a single contribution 
oc ( m ^d ■ Moreover, since this case is exactly soluble, we have access to all orders of the derivative expansion. Indeed, 
for the inhomogeneous magnetic background (38), the inverse mass expansion (28) becomes a double sum 



8tt 3 / 2 



E 



(mAwE 

v > k=2 



2eB 

nrrt 2 



2 k 



r(2fc + j)r(2fc + j-2)g 2fc+2j 

j\{2k)T{2k + j + \) 



(39) 



with the two expansion parameters being the derivative expansion parameter, and the perturbative expansion 
parameter, ^|. Note that the expansion coefficients are known to all orders, and are relatively simple numbers, 
just involving the Bernoulli numbers and factorial factors. It has been checked in [38] that the first few terms of 
this derivative expansion are in agreement with explicit field theoretic calculations, specialized to this particular 
background. 

Given the explicit series representation in (39), we can check that the series is divergent, but Borel summable. 
This can be done in several ways. One can either fix the order k of the perturbative expansion in (39) and show 
that the remaining sum is Borel summable, or one can fix the order j of the derivative expansion in (39) and show 
that the remaining sum is Borel summable. Or, one can sum explicitly the k sum, for each j, as an integral of a 
hypergeometric function, and show that for various values of the remaining derivative expansion is divergent 
but Borel summable. These arguments do not prove rigorously that the double series is Borel summable, but give a 
strong numerical indication that this is the case. It is interesting to note that an analogous double-sum structure also 
appears in the renormalon-OPE analysis in the talk by E. Gardi in these Proceedings. 

It is also possible to compute the closed-form effective action for a time dependent electric field pointing in the 
z-dircction: 



E(t) =Escch 2 (t/ T ) 



(40) 



where r characterizes the temporal inhomogeneity scale. A short-cut to the answer is to note that we can simply make 
the replacements, B 2 — » — E 2 , and A 2 — > — r 2 , in the magnetic case result (39). In particular this has the consequence 
that the alternating divergent series of the magnetic case becomes a non-alternating divergent series, just as was found 
in the Euler-Heisenberg constant-field case. Fixing the order j of the derivative expansion, the expansion coefficients 
behave for large k (with j fixed) as 



„0") 



(-iy+ k T(2k + 3)T(2k + 3 + 2)B 2k+2]+2 T(2k + 3j - ±) 



r(2fc + 3)r(2fc + j + |) 



(27r) 2 J'+ 2fc + 2 



Note that these coefficients are non-alternating and grow factorially with 2k, as in the form of (24). 
Borel dispersion formula (26) gives 
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Remarkably, this form can be resummed in j, yielding a leading exponential 



m 
8V 3 



ImS ^ — — « I — o 



eE 
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III 



5/2 



exp 



m 2 7r J 
~~e~E \ 



^ 1 / m \ 2 
~ 4 \~eEh 



(41) 



Applying the 



(42) 



(43) 



We recognize the first term in the exponent as the leading exponent in the constant field result (32). Thus, the second 
term may be viewed as the leading exponential correction to the constant-field answer (32). This is exactly what 
we set out to find, and we see that it arose through the divergence of the derivative expansion. I stress that this 
exponential correction is not accessible from low orders of the derivative expansion. This gives a Dyson-like argument 
that the derivative expansion must be divergent, since if it were not divergent, there would be no essential difference 
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between the electric and magnetic cases, and there would be no correction to the exponent of the imaginary part of 
the effective action. However, we know, for example from WKB, that there is such a correction, and so the derivative 
expansion must be divergent. 

In fact, the situation is even more interesting than the result (43) suggests. We could instead have considered doing 
the Borel resummation for the j summations, at each fixed k. Then for large j, the coefficients go as 

jfc) _ ( i y+fc r(j + 2fc)r(j + 2k - 2)B 2k+23 9_ 2fc r(2j + 4fc-|) 

j { ' T{j + l)T{j + 2k+\) ~ (27r)2j+2fc I ) 

which once again are non-alternating and factorially growing. Applying the Borel dispersion formula (26) gives 

4^/2 ( 2 jfe)! 6 (45) 
Once again, this leading form can be resummed, yielding 

eEr 



m 3/2 

^~^-^exp 



-2txvht [ 1 



(46) 



Note that this leading exponential form of the imaginary part is different from that obtained in (43), and moreover, 
it is different from the constant-field case (32). The resolution of this puzzle is that there are two competing leading 
exponential behaviours buried in the double sum (39), and the question of which one dominates depends crucially on 
the relative magnitudes of the two expansion parameters, the derivative expansion parameter, and the perturbative 
expansion parameter, Another important quantity is their ratio, since this sets the scale of the corresponding 
gauge field: 

A(t) e£>rtanh(i/T) eEr eE/m 2 , Ar7 . 

= ^ = (47) 

m m m l/(mr) 

Thus, it is natural to define a "nonperturbative" regime, in which 3> 1. Then mr » 2^, so that the dominant 

exponential factor is exp[— ^] 3> exp[— 27rmr]. In this regime, the leading imaginary contribution to the effective 
action is given by the expression (43), and we note that it is indeed nonperturbative in form, and the correction in 
the exponent is in terms of the small parameter <C 1. On the other hand, in the "perturbative" regime, where 

<C 1, this means that tot <C f^, so that the dominant exponential factor is exp[— 2-kvht\ » exp[— In 
this regime, the leading imaginary contribution to the effective action is given by the expression (46), and is in fact 
perturbative in nature, despite its exponential form. 

These results are completely consistent with the WKB approach developed by Brezin and Itzykson [34] and Popov 
[35]. For a time-dependent electric background E{t) = A z (t), in the z-direction, the WKB expression for the imaginary 
part of the effective action is: 



ImS ~ J d 3 k exp[-7rO] 



(48) 



where = ^ / ^Jm 2 + k\ + (k z - eA z (t)) 2 . Applying this WKB analysis to the (exactly soluble) case E(t) = 

£sech 2 (t/r), one obtains precisely the leading results (43) or (46), depending on whether we are in the "non- 
perturbative" » 1, or "perturbative" <C 1 regime. This serves as a useful cross-check of the somewhat 
formal Borel analysis. 



III. PERTURBATIVE - NONPERTURBATIVE DUALITY IN QES SYSTEMS 

In this last section I discuss some recent results [39] concerning a new type of perturbative - nonperturbative 
connection that has been found in certain special quantum mechanical systems. We benefited from discussing these 
systems with Arkady, and I hope he enjoys the results! 

Quasi-exactly solvable (QES) systems are those for which some finite portion of the energy spectrum can be found 
exactly using algebraic means [40]. A positive integer parameter J characterizes the 'size' of this exact portion of the 
spectrum. Two simple examples are : V = x 6 — (4J — l)x 2 , and V — sinh 2 x — (2 J — 1) coshrr. For a QES system 
it is possible to define a quadractic form, H — J2 a b c abJ a Jb + J2 a da J a , in terms of si (2) generators of spin J, such 
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that the eigenvalues of the algebraic matrix H are the QES eigenvalues of the original system. It is interesting that 
algebraic hamiltonians of this form are widely used in the study of tunneling phenomena in single-molecule magnets 
[41]. 

In [42], the large J limit of QES systems was identified as a semiclassical limit useful for studying the top of 
the quasi-exact spectrum. It was found that remarkable factorizations reduce the semiclassical calculation to simple 
integrals, leading to a straightforward asymptotic series representation for the highest QES energy eigenvalue. The 
notion of energy- reflection (ER) symmetry was introduced and analyzed in [43]: for certain QES systems the QES 
portion of the spectrum is symmetric under the energy reflection E <-> — E. This means that for a system with 
ER symmetry, there is a precise connection between the top of the QES spectrum and the bottom of the spectrum. 
Coupled with the semiclassical large J limit, the ER symmetry therefore relates semiclassical (nonperturbative) 
methods with perturbative methods [43] . In this section I discuss a class of periodic QES potentials for which the ER 
symmetry is in fact the fixed point (self-dual point) of a more general duality transformation. The duality between 
weak coupling and semiclassical expansions applies not just to the asymptotic series for the locations of the bands 
and gaps, but also to the exponentially small widths of bands and gaps. 

Consider the quasi-exactly solvable (QES) Lame equation [44]: 

+ J(J + l)vsn2 ^ v) - \ J{J + 1} } *w = E * ^ • (49) 

Here sxi(4>\v) is the doubly-periodic Jacobi elliptic function [27, 44], the coordinate <fi £ R 1 , and E denotes the energy 
eigenvalue. The real elliptic parameter v lies in the range < v < 1. The potential in (49) has period 2K(y), where 
K{y) is the elliptic quarter period. The parameter v controls the period of the potential, as well as its strength: sec 
Fig. 7. As v — > 1, the period 2K{v) diverges logarithmically, 2K{v) ~ In (y^), while as v — > 0, the period tends to a 
nonzero constant: 2K(y) — ► w. In the Lame equation (49), the parameter J is a positive integer (for non-integer J, 
the problem is not QES). This parameter J controls the depth of the wells of the potential; the constant subtraction 
— \ J(J + 1) will become clear below. 




FIG. 7: The potential energy in (49) as a function of (f>, for 3 = 2. The solid curve has elliptic parameter v = 0.95, for which 
the period is 2if(0.95) « 5.82. The dashed curve has v = 0.05, for which the period is 27^(0.05) fa 3.18. Note how different 
the two potentials are; and yet, their spectra are related by the duality transformation (53). 

It is a classic result that the Lame equation (49) has bounded solutions *&(4>) with an energy spectrum consisting of 
exactly J bands, plus a continuum band [44]. It is the simplest example of a "finite-gap" model, there being just a finite 
number, J, of "gaps" in the spectrum. This should be contrasted with the fact that a generic periodic potential has an 
infinite sequence of gaps in its spectrum [45]. We label the band edge energies by E\, with I — 1,2,..., (2J+ 1). Thus, 
the energy regions, E21-1 < E < E21, and E > E21+1, are the allowed bands, while the regions, E21 < E < E21+1, and 
E < E\, are the gaps. 

Another important classic result [46, 47] concerning the Lame model (49) is that the band edge energies Ei, for 
I = 1, . . . , 2 J + 1, are simply the eigenvalues of the finite dimensional (2J + 1) x (2J + 1) matrix 

H = Jl + vJl- l -J{J + 1)1 (50) 

where J x and J y are su(2) generators in a spin J representation and I is the unit matrix. Thus the Lame band edge 
spectrum is algebraic, requiring only the finding of the eigenvalues of the finite dimensional matrix H in (50). For 
example, for J = 1 and J = 2, the eigenvalues of H are: 



J = 1 : £1 = -l + u, 

E 2 = 0, 
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(51) 



J = 2 



E2 
E3 
E A 

E 5 



-1 + 2v-2y/l-v + v 2 , 
-2 + i/, 
-2 + 4!/, 
1 + z/, 



-1 + 2^ + 2 VI - f + v 2 



(52) 




0.2 



0.4 



0.6 



FIG. 8: The energy bands (52) for the J = 2 Lame system (49), as a function of the elliptic parameter v. The shaded areas on 
this plot are the bands, while the unshaded areas are the gaps. The top band actually continues up to E — > 00. 



The spectrum of the Lame system (49) has a special duality: 

E[v] = -E[l - 



(53) 



That is, the spectrum of the Lame system (49), with elliptic parameter is, is the energy reflection of the spectrum of 
the Lame system with the dual elliptic parameter \ — v. In particular, for the band edge energies, Ei, which are the 
eigenvalues of the finite dimensional matrix H in (50), this means 



EiM 



-E. 



2J+2- 



l = 1,2,..., 2.7+1 



(54) 



This duality can be seen directly in the eigenvalues of the J — 1 and J = 2 examples in (51) and (52). The proof of 
the duality (54) for the band edge energies is a trivial consequence of the algebraic realization (50), since 
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(Jl + (1- V )jl)--J{j+I)l 



(55) 



Noting that [J| + (1 — v) J 2 ] has the same eigenvalues as [J 2 + (1 — v)J 2 ] , the duality result (54) follows. It is instructive 
to see this duality in graphical form, in Fig. 8, which shows the band spectra as a function of v. The transformation 
v — > 1 — v, with the energy reflection E — > — E, interchanges the shaded regions (bands) with the unshaded regions 
(gaps). The fixed point, v = ^, is the "self-dual" point, where the system maps onto itself; here the energy spectrum 
has an exact energy reflection (ER) symmetry. 

In fact, the duality relation (53) applies to the entire spectrum, not just the band edges (54). This is a consequence 
of Jacobi's imaginary transformation [27]. Making the coordinate transformation 



i{<j> — K — i K') , 



the Lame equation (49) transforms into 

d 2 



f/2 



+ J(J + 1) (1 - v) sn 2 (0'|l -v)- -J(J + 1) } *(<£') 



(56) 

-E^{4>') 

(57) 

So solutions of (49) are mapped to solutions of the dual equation (57), with v — > 1 — v, and with a sign reflected 
energy eigenvalue: E — > —E. 

To see why bands and gaps are interchanged under our duality transformation (56), recall that the independent 
solutions of the original Lame equation (49) can be written as products of theta functions [44], and under the change 
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-4-2024 -4-2024 



FIG. 9: Contour plots in the <f> plane of the real (left) and imaginary (right) parts of the potential function sn {<t>\v), for 
v = |. Note that the potential is real along the real <fi axis, and along the line Re(4>) = K(^) « 1.85 (as well as the periodic 
displacements) . This corresponds to the rotation (56) . 



of variables (56), these theta functions map into the same theta functions, but with dual elliptic parameter. However, 
they map from bounded to unbounded solutions (and vice versa), because of the "i" factor appearing in (56). Thus, 
the bands and gaps become interchanged. 

As an aside, I mention that the Lame system plays a distinguished role in the theory of su(2) BPS monopoles. For 
example, Ward has shown [47] that the Lame equation factorizes, using the Nahm equations (which are fundamental 
to the construction of monopole solutions) 

^-=T^e jk i[T k ,T l ] (58) 

To see this, define the quaternionic operators 

A± = l 2n ^- s ±T 3 {s)®a J (59) 
where Tj(s), for j = 1, 2, 3, take values in an n-dimensional represenation of the Lie algebra su(2). Then 

A±A T = (l n -^ - TjT^J ® 1 2 + ( T § - \e jkl [T k , T,]) ® a, (60) 

and this combination is real if Tj satisfy the Nahm equations (58). Furthermore, using the solution T\ = — Y / I'sn(s) t\, 
T 2 = i^/v ca{s)t2 1 T3 = zdn(s)t3, where tj are su(2) generators, it follows from properties of the elliptic functions 
that 

TjTj = (t? + 4 + t 2 3 )vsn 2 ( S ) -{vt 2 2 + t§) (61) 

Thus, (60) provides a factorization of the Lame operator. Subsequently, Sutcliffe showed [48] that the spectral curve 
for the Nahm data for a charge n = 2j + 1 su(2) monopole is related to a j-gap Lame operator, and corresponds 
physically to 2j + 1 monopoles aligned along an axis. 

Returning to the perturbative-nonperturbative duality (53), I first discuss how this operates for the locations of 
the bands and gaps in the Lame spectrum. The location of a low-lying band can be calculated using perturbative 
methods, while the location of a high-lying gap can be calculated using semiclassical methods. The exact duality 
(53) between the top and bottom of the spectrum provides an explicit mapping between these two sectors. Defining 
k = yf J {J + 1), we see that 1/k is the weak coupling constant of the perturbative expansion. Simultaneously, 1/k 
plays the role of Ti in the quasiclassical expansion. 

In the limit J — > 00, the width of the lowest band becomes very narrow, so it makes sense to estimate the 
"location" of the band. In fact, the width shrinks exponentially fast, so we can estimate the location of the band to 
within exponential accuracy using elementary perturbation theory. A straightforward calculation [39] shows that the 
lowest energy level is 



k 2 k 2 8V^« 3 W 



(62) 



We now consider the semiclassical evaluation of the location of the highest gap, in the limit J — > 00. First, note 
that for a given v, as J — > 00 the highest gap lies above the top of the potential. Thus, the turning points lie off the 
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real (f> axis. For a periodic potential the gap edges occur when the discriminant [45] takes values ±1. By WKB, the 
discriminant is 



A(E) = cos 



S n (P) 



(63) 



n=0 



where P is the period, and S n {x) are the standard WKB functions which can be generated to any order by a simple 
recursion formula [16]. The J th gap occurs when the argument of the cosine in the discriminant (63) is Jn: 



°° 1 / 1 1 1 \ 

kV — S n (2K) = Jn = ttk[1- — + —r - — -. + ...) 

^ K n v ' \ 2k 8k 2 128k 4 / 



(64) 



where 2K is the period of the Lame potential, and where on the right-hand side we have expressed J in terms of the 
effective semiclassical expansion parameter 1/n. This relation (64) can be used to find an expansion for the energy of 
the J th gap by expanding 



oo 



-2 ' 



(65) 



The expansion coefficients £; are fixed by identifying terms on both sides of the expansion in (64). A straightforward 
calculation [39] leads to 



1 , , 2-v (-2 + 2v + v 2 ) 

E = — K — K + — — + ± —== ' 

2 4 16kVT^17 



(66) 



Comparing with the perturbative expansion (62) we see that the semiclassical expansion (66) is indeed the dual of 
the perturbative expansion (62), under the duality transformation v — > 1 — v and E — > — E. 

This illustrates the perturbative - nonperturbative duality for the locations of bands and gaps. However, it is even 
more interesting to consider this duality for the widths of bands and gaps, because the calculations of widths arc 
sensitive to exponentially small contributions which are neglected in the calculations of the locations. 

The width of a low-lying band can be computed in a number of ways. First, since the band edge energies are 
given by the eigenvalues of the finite dimensional matrix H in (50), the most direct way to evaluate the width of the 
lowest-lying band is to take the difference of the two smallest eigenvalues of H. This leads [49] to the exact leading 
behavior, in the limit v — > 1, of the width of the lowest band, for any J : 



A ^algebraic 
Aft band 



8JT(J + 1/2) 

4^r(j) 



J -I 



(67) 



This clearly shows the exponentially narrow character of the lowest band. 

In the instanton approximation, tunneling is suppressed because the barrier height is much greater than the ground 
state energy of any given isolated "atomic" well. The instanton calculation for the Lame potential can be done in 
closed form [49], leading in the large J and v — * 1 limit to 



A z^instanton 
A ^band 



8J 3 / 2 
0F4 J 



(68) 



which agrees perfectly with the large J limit of the exact algebraic result (67). Thus, this example gives an analytic 
confirmation that the instanton approximation gives the correct leading large J behavior of the width of the lowest 
band, as v — > 1. 

Having computed the width of the lowest band by several different techniques, both exact and nonperturbative, we 
now turn to a perturbative evaluation of the width of the highest gap. First, taking the difference of the two largest 
eigenvalues of the finite dimensional matrix H in (50), it is straightforward to show that as v — > 0, for any J, this 
difference gives 



AE, 



algebraic 
gap 



8jr(j + 1/2) 
4^r(j) 



i 



j -i 



■v + . 



(69) 



which is the same as the algebraic expression (67) for the width of the lowest band, with the duality replacement 
v — > 1 — v. But it is more interesting to try to find this result from perturbation theory. From (69) we see that 
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the width of the highest gap is of J th order in perturbation theory. So, to compare with the semi-classical (large J) 
results for the width of the lowest band, we see that we will have to be able to go to very high orders in perturbation 
theory. This provides a novel, and very direct, illustration of the connection between nonperturbative physics and 
high orders of perturbation theory. 

It is generally very difficult to go to high orders in perturbation theory, even in quantum mechanics. For the Lame 
system (49) we can exploit the algebraic relation to the finite-dimensional spectral problem (50). However, since H 
in (50) is a (2J + 1) x (2J + 1) matrix, the large J limit is still non-trivial. Nevertheless, the high degree of symmetry 
in the Lame system means that the perturbative calculation can be done to arbitrarily high order [39]. The result is 
that the splitting between the two highest eigenvalues arises at the J th order in perturbation theory, and is given by 

A pport. theory _ 8 ( 2 -0 ! J _ 8JT(J + 1/2) j 

-457 pTIjif " - 4V^r(j) " (70) 

This is in complete agreement with (69), and by duality agrees also with the nonperturbative results for the width 
(67) of the lowest band. 



IV. CONCLUSIONS 



To conclude, there are many examples in physics where there are divergences in perturbation theory which can be 
associated with a potential instability of the system, thereby providing an explicit bridge between the nonperturbative 
and perturbative regimes. While this is not the only source of divergence, it is an important one which involves much 
fascinating physics. Arkady has made many important advances in this subject. On this occasion it is especially 
appropriate to give him the last word: 

"The majority of nontrivial theories are seemingly unstable at some 
phase of the coupling constant, which leads to the asymptotic 
nature of the perturbative series." A. Vainshtcin, 1964 [1] 
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